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I. INTRODUCTION 

InLegro-differential equations arise in many engineering and scientific disciplines, 
often as approximations to partial differentia! equations, which represent much of 
the continuum phenomena. Many forms of these equations are possible. In this re- 
port, we study identification of integro- differential equations, proposed to describe 
unsteady aerodynamics and aeroelasticity phenomena. The forms of equations stud- 
ies in this report are closely related to the research of Tobak and Schiff [ 1-4 J among 
others. 

Let us consider the symmetrical longitudinal dynamics of an aircraft. In the 
linear quasi-steady regime, the aerodynamic pitch-moment coefficient, C,„ , can be 
written as a linear combination of an gle-of- attack, a , pitch rate, q , and elevator 
deflection, S e . 


C„i — C 1 ,,,., + O mo a + C,„ t q + C mji 6 c (l) 

The equation for the pitch coefficient implies that its values at any point in time 
depends only on the instantaneous values of a , q , and S e . This is a good approx- 
imation in quasi-steady flow. In regions where the flow field is highly unsteady, or 
if the elevator deflections are changed rapidly, it is reasonable to believe that the 
pitch moment coefficient depends on the past as well as the instantaneous values 
of the flow variables. Such conditions may occur in separated flow regions such as 
high angle-of- attack, stall, or spin. A reasonable extension to include past history 
of o in Equation (l) would lead to: 

r 

C, n — Cm,, + O m + J (r)c*(£ — r)dr + C m ^ q -f- C m „ t (2) 

o 

Of course, q and S r terms could be expanded similarly. Using the above equation 
for pitch movement coefficient, the pitch rate dynamics would be as follows: 

T 

fq - q^sc{C, n „ + C mn a + J F,, o (r)a(t - r)dr + — C m< q C m , e <5 e } (3) 

o 

where : 

I : Pitch moment of inertia 

q : Dynamic Pressure 

s : Reference area 

c : Reference wing chord 

V : Aircraft Speed 

We, thus, have the simplest form of integro-differential equation arising in un- 
steady aerodynamics. Similar results have been obtained by more systematic aero- 
dynamic analyses. The integral term may also be written in terms of d as it is 


usually done in aerodynamic analysis. Suppose we write ( assuming a = 0 for 

t - oo ) .* 


C,r. - C„ lo + J /I,„ A (r)a(t - t)Ut 4- + C, nt> 6 e (-1) 

0 

The form has been referred to as the indicial representation by aerodynamicists [3j, 
This equation is equivalent to Equation (2). Integrating the a term by parts : 


oo CO 

J H<h A (r)a(i - r)dr ~ -//,„ A ( r)a(t - r)j™ 4- J , ( I - r)a(t - r)dr 

0 0 

co 

= II, (O)or(i) 4- J - r)a(t - r)dr. 


(5) 


Thus, Equation (4) is equivalent to Equation (3), if 

= C mo 

and 


~II mA - K>„ 0 < r < T 

= 0 T < r (fl) 

ft is reasonable to assume that in quasi-steady flow /f,„ d ( r ) will reach a constant 
value as r -* oo . In practice, it is often convenient to use one form for identi- 
fication and convert it into the other form for analysis and understanding of the 
aerodynamics phenomenon involved. 

In unsteady aerodynamics and aeroelasticity, Equation (3) could be generalized 
in any of the following ways and possibly others. 

( a ) C m ,, , and C„ ti ' could be nonlinear functions of angle-of-attack, pitch 
rate and elevator deflection. 

(b) The equation for C m could involve double or higher order integrations, in- 
volving products of independent variables. 

(c) Full six degree-of-freedom dynamics would bring in nonlinear kinematic 
terms. 

This report addreses the following problems: 

(a) Section 2 develops techniques for identifying integro- differential equations 
models from test or simulation data. 

(b) Section 3 shows when the the integro-differential equation may be approxi- 
mated by an ordinary differential equation. 

(c) Conditions under which the integral term may be identified accurately are 
analyzed in Section 4. 

The concepts developed in the previous sections will be applied to simulation 
data and high angle-of-attack test data in Section 5. Section 6 gives summary and 
conclusions of the results presented in this report. 
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II. IDENTIFICATION OF INTEGRO-DIFFERENTIAL EQUATIONS 


The problem of estimating aircraft stability and control parameters has been 
researched extensively and many solutions have been proposed [6-8]. The identifi- 
cation of intcgro-difTerential equation models poses difficulties such that previous 
methods need to be extended. The discussion in this section will be limited to the 
form of Equation As has been shown in the previous section, Equation (3) is 
equivalent to tin ;rm involving a in the integral. Equation (3) is more suitable for 
identification, because a is usually measured while a is not. 

The identification problem for the pitch moment equation (3) involves the de- 
termination of (a) parameters C mo ,C ina ,C m ,and C,„ it , (b) function F m<i (r) 0 < r < T, 
and (c) T. For the purpose of our discussion, we will assume that measurements 
<7(0, Q (0, andfi e (0 are available from a flight test at sample points t k ,k - 1,2.. .A. Other 
measurements can be included in the identification process without changing the 
procedure significantly. 

The problem of identifying integro-differential equation models is, thus, different 
from the parameter estimation problem because it is necessary to identify a function 
in addition to parameters. There are two approaches to address the identification 
problem - identify the function directly, or approximate the function. In this section, 
we discuss each of the approaches. 

2.1 D irect Identification of Integro-Differential Equations 

To estimate the F„ la (t) and the parameters, we consider Equation (3) and min- 
imize the following performance index by choosing C,„ 0 ,C ma ,C mii , and as well 

35 F„,„ (r) . 


J = \ J2 (<?">(*) - 

fc=i 

The first gradient of J with respect to F„ ln (r) is computed as follows: 


(7) 


9F{t) 


dJ 
d F,n„ 


N 


" “ H ~ g{k)) 


fc=i 


3q 

BF m . 


and 


d r 3q 

dt { dF mo 


} = —r-Wt- 


t) A + 


3q 

9F m . 


( 8 ) 


(9) 


where A represents a small region around r over which F ma (r ) is perturbed. Note 
that because F m „(r) is a continuous function, perturbation in F m<1 (r) at a single 
point will not change the value of q at all, since that will not change the integral. 
To perform any unique identification at all, it is necessary to assume that F jna {r) 
possesses certain continuity properties. This is in any case a good assumption for 
most physical systems. 


Since a ancl q are highly nonlinear functions of F,„ a (r) , an iterative procedure 
is needed to optimize the likelihood function. Either a first or second gradient 
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procedure can be used. In a first gradient approach, we could start from a initial 
value of F nto (r) and update it as follows ( k. represents iteration number) 


- F„, n (r)U £ 3r(T)|fe (10) 

The value of £ is selected based on a line search. The procedure can be repeated 
many times until convergence occurs. A second gradient approach could also be 
used. 

The determination of the first gradient can be computationally time- consuming. 
Note that the first gradient is a function of r and must be evaluated for all values 
of r. That would require differential Equation (9) to be solved for all r. Second 
gradient would be even more difficult to compute. Thus, it would be desirable from 
a computational point of view to divide up the interval between 0 and T into M 
intervals and estimate the values of (r) at M + 1 individual points, e.g,, 

fmji—), « = 0 , 1,2 ( 11 ) 

This would require the solution to M + 1 equation of the form of (6) to get the 
first gradient. To cut down the computation time, it is also possible to start with 
fewer intervals in early iteration and increase the number of intervals as the problem 
converges. 

The problem of identifying T can be handled in one of two ways. One approach 
is to select a large enough T and estimate F m „(r) . In principle, the estimated values 
of F„ ln (r) for r greater than the true value of T should be small or zero. Else an 
estimate could be obtained by computing the gradients of the cost functional as 
follows: 


8J f, ,,u 5 9( fc ) 

9T = dT = ■■■ 

L J k=l 


dT 


( 12 ) 


£r^i\ _ 

dear’ 




{F, no (T)a(t 


T) + C m , J^} 


( 13 ) 


The joint estimation of C7 m „, C mo , C m> , and C r „ 4t , F„ la (r) , and T is a straight 
forward extension of the above procedure. If ^(i-j is estimated at many points, a 
first order gradient procedure may be the only viable option. For estimating a few 
points on the /' , mo (r) profile, a second gradient procedure is also usable. 


2.2 Identification of Integro-Differential Equation Models by Approximation 


Assumptions for Approximation 

The identification problem is considerably simplified if the integral term is ap- 
proximated by a sum of predefined functions, called basis functions. These approxi- 
mations convert unknown functions into a small number of parameters, but require 
the following assumptions. 

1. F ma (r) aid other unknown functions have small high frequency components 
(note thi t since a multiplies F m(j (r) , the integral term in Equation (3), can 
be a major contributor to the pitching moment at high frequency). 


2. F m „(r) an d other unknown functions have finite tails; i.e., are zero for r > T. 

3. F m „(r) is a continuous function, which converges quickly with any reasonable 
set of basis functions. 

This identification procedure, therefore, consists of two steps - representation of 
indicia! function T rriii (r) by a set of basic functions and associated parameteis, and 
then estimation of the resulting parametric models. 

Representation of Integro- Differential Terms 

Smooth functions have been parameterized in many different ways. For the 
purpose of identification, the following forms appear most appropriate. 

(l) Polynomials : F,„ a (r) is written as a polynomial in r 


F,n„ (r) pa At) -I- A I t -I- A 2 t 7 + 

« 0 


0 < r < T 
T < r 


( 14 ) 


Ad, A lt A 2l A a are the unknown parameters. Numerical conditioning can be 
improved by either 

(a) expanding the polynomial in , or 

(b) using orthogonal polynomials. 

(2) Splines : Linear, quadratic on cubic splines are excellent approximation 
to most continuous functions. Splines are polynomials whose coefficients 
change at certain break points, called knots [9]. The coefficients change in 
a constrained way to ensure continuity at knots. The unknown parameters 
are the polynomial coefficients as well as knot locations. 

B-splines are preferred from a numerical viewpoint [9] . 

(3) Impulse Response of Rational Models: This form is useful when T — * oo, 
since any other form requires too many parameters. The Laplace transform 
of F,„ a (rj could be written as a rational function of s 


*n» 


m*) 

D( s ) 


( 15 ) 


To ensure that F mti (r) -* 0 as r — * oo, £>(s) must have stable poles. For 
example ) T, t i „ ( t ) might be written as : 



( 16 ) 


t 

J F m „(t - r)a(r)dr w + c 2 f 2 

— OO 


( 17 ) 


This will correspond to 


N{») = + b 2 

D(s) s + ai 3 + a 2 

This form has been used in aeroelastic studies. 


( 18 ) 
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(<l) Slops : The indicial functions could be wiiu i as 


= “1, 

0 < r < - 

- 

n 


7 

2 7 

= 0 2 , 

— < T < 

n 

n 

= a r»p 

(n - i)7 

< T <7 


n 


= 0, 

T <T 



(19) 


Identification Procedure 

Once the integro-difTercntial equation model is converted into a parametric form, 
previous estimation procedures can be used to determine the parameters as well as 
the model form [6]. The entire procedure is shown schematically in Figure 1. 
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HI. APPROXIMATION OF INTEGRO - DIFFERENTIAL EQUATIONS 


In this section, errors in approximating intcgro-difFerential equations by ordi- 
nary differential equations are estimated. The analysis addresses the specific class 
of equations encountered in unsteady aerodynamics and aeroclasticity. 

Consider the integral term in Equation (2). Under certain conditions, this term 
is approximated as an add on to C, na . In the current analysis, we assume that the 
approximated form represents the entire pitch momert coefficient due to angle-of- 
attack. 


r 

f F m Jr) a (t - T )dT = C„ >a a(t) (20) 

0 

where C m „ is a const. \nt . For the two sides to equal in steady-state 

T 

f F'nJr)dT = C m „ (21) 

0 

The approximation error is 

r 

e - j F m „( r )a (t - r)dr - C m „a(t) (22) 

o 


In the discussion of the approximation, we assume that a(i) is a continuous 
function with no jumps, such that a(t) can be defined at each point. Since the 
secant slope cannot exceed the maximum tangential slope in an interval, we get : 

||q(£ - r) - a(t)|| < r||d m , 1I (£)|| (23) 

where a„ lax (t) is the maximum value of a in the time interval from i.-t to ( . Equation 
(22) therefore can be written as : 


r 

£ = J F m „(r)(a(t - r) - c*(£))cir 

o 

T 

< j ||F mo (r)!|( Q (£-r)-a(t))|ldr 
0 

T 

< <im„x(t) J ||F m „(r)Hrrfr (24) 

0 

Thus, the maximum error depends on a max (t) and the first moment of F m „(r) . 
If we make no assumption about the nature of F mo (r) , then 
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< < Ta, ltai {l)C m 


General 


(2fi) 


I'or rectangular and triangular F„ la (r) , the approximation takes the forms 

e < ^Ta max (t)C ma : rectangular 

^ ^TFlHI triangular 

Since the maximum approximated value of the integral term is C„,„ (t), the 

relative error becomes 


< T- 


d!) 

.(0 


< ^ fp Griuix (t) 

“ 2 a„ 


< 

3 a 


max 

&max 


(t) 

{<-) 


.(0 


General C m „(r) 
rectangular 
triangular 


( 20 ) 


A similar result is obtained in approximating Equation (4) (Appendix A) 

I hese error equations can be further simplified for improved understanding for 
motion at a single frequency w. In that case 


Qmm ^ 

So the relative error may be written as : 


(27) 


c r <u>T : General (r) 

< -wT 1 • rectangular 

< ~wT : triangular 

O 

The error introduced by the approximation increases with frequency of the motion 
and the maximum time delay for which the integral term is significant. 


What do these two terms depend on? The frequency of the motion varies with 
applied inputs, natural short period dynamics and external disturbances. Most 
fixed-wing aircraft appear to have the natural frequency of the short period dy- 
namics in the neighborhood of 0.5 Hz, though this number could be diiTerent for 
nonconventional aircraft. Pilot applied inputs are typically less than 1 Hz. Most 
turbulence spectra is also below 1 Hz. 

Two factors can increase w significantly. First, unsteady aerodynamics, vortex 
shedding and stall/spin can cause large rapid motions. Secondly, the aeroelastic 
phenomenon, because of high structural frequencies, almost always occurs at higher 
values of w. 


The Lime T can be physically thought of as the delay in reestablishing a new 
flow field following a change in one or more aerodynamic variables. Under quasi- 
steady conditions, T is proportional to the time taken by the flow field to travel a 
characteristic distance L. Thus 


T L y ( 28 ) 

where V is aircraft speed. Thus T increases at low aircraft speed. If there is little 
or no interaction between the wing and the La!', the characteristic distance is the 
wirg chord, while it could be as large as the aircraft length if vortex shedding at 
the nose can impact the tail. 

bet us evaluate three ifitualions to determine if the integral term can be impor- 
tant in evaluating pitching moment coefficient (in each case assuming the integral 
term is no v/orse than a triangular form). 

Case ■ A fighter aircraft 60’ long 10’ wing chord, traveling at 800 ft /sec . Short 
period natural frequency is 0.5 Hz 


800 80 

c r < - * 2 * «• * .S * — — .013 
~ 3 80 

Case 2: The same aircraft in stall, nose/wing vortices hitting tail, forward speed 120 
ft sec , forced dynamics at 0.75 Hz 

1 00 
c r < - * 2 * n * .7G * — - = .78 

“ 3 120 

Case 3: The same aircraft with flutter behavior near 4 Hz. Forward speed 1000 ft 
sec . Unsteady aerodynamics requiring twice the normal time to settle down 

1 10 

c r < -*2*ir*d* 2-—- = .17 

3 1000 

The integral terms could be important in the last two cases. It has been conjectured 
by Tobak and SchifT [4] that in the quasi-steady state region at low anglc-of-attack, 
the indicial terms are not important. 

It should be noted that applying high frequency inputs does not necessarily 
make the indicial term more important. At high frequency, large amplitude inputs 
are needed to produce small motions in a and q . Thus, even though the indicial 
term is significant compared to C ma a, both of these terms could be small compared 
to the term. 

The main result of this section is therefore, as follows. 

Tl it integro- differential equations arising in unsteady aerodynamics and aeroelasticity may be 
approximated by differential equations when 


wT « 1 

where w is the natural frequency of the motion and T is the time delay over which the indicia/ 
function is significant. 
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IV. IDENTIFIABILITY 

To study the identifiability problem, wo will start with a more general model 
for the pitch moment equation 


CO 

-Qoa a c { £«>„ -|- C rI1<1 q I q -|- C„, a . &, + J — r)dr 

o 


oo <x> 

' 2V / T ) dr + J F mtt 6 e [l - r)dr) 

o » 


( 29 ) 


The identifiability of the single-axis planar motion is analyzed to indicate some 
basic limitations to obtaining the values of functionals in integro-differcntial models. 
Let us assume that the a(<) and <5,(0 can change independently. Correlation between 
»(<) and , caused by aircraft dynamics, will further reduce identifiability. It is 
also assumed that reference values of , and a(£) and 6 r (t) are selected such the 
„ is zero. 

Because of superposition, the identifiability of linear systems can be determined 
by considering a set of inputs at various frequencies and amplitudes. Let us write 
Equation (24) in the frequency domain 


jwq(]u} - + c m<l a{]w) -f- ^yC m ,q[jw) + C mi '6 e ( jut) 

-< F >"n (jw)or(yw) + (joj)q(ju) -»• F, ni ' (jw )5 e (jw)) (30) 

We will apply sinusoidal inputs to the aircraft . Let 

a = a 0 e J “' 

6n = SoS* 

To evaluate identifiability, we can select two linearly independent combinations 
of a„ and S„ , e.g., 


a„ — 1 , S 0 — 0; and 

a„ = 0 , S„ = 1 . ( 31 ) 

Because of linearity, other values of a 0 and S 0 do not afTect system identifiability even 
though they lead to reduced estimation errors in parameters which are identifiable. 
Thus at any w there are four scalar equations (two complex equations) and nine 
unknowns ((?,„,, C ma ,C mit and real and imaginary parts of F ma , F„ u and F mit . The 
parameters C m „, £7,,,,, C t „ 4t are not functions of frequency. 

At zero frequency, we have two equations and three unknowns ( C m „, C m f , C,„ 4< ) . 
The zero frequency input data provides a mechanism to solve for C m<j and <? m ,in 
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terms of ( or any two variables in terms of the third). Thus, at any frequency 
there arc four equations and seven unknowns. Without further assumptions, linear 
indicial function representation in the pitch plane !b not completely identifiable from 
flight test data. 

The inclusion of C’„, u in terms of C,„ 4< by including a Lest point at 0 and 
A„ --- o in addition to the two test points described by Equation (20). 

The non-identifiability of F m ft) , (t) , and F„ u> (t) indicates that further con- 

straints arc needed on the general nature of these functions. Some of the properties 
which should hold because of the basic naLure of fluid flow can be used to bring 
about idenlifiability. These properties relate to smoothness condition discussed in 
Section 2. 




V. SIMULATION AND FLIGHT TEST RESULTS 

Many inlegro-differential models ha’ e been identified using the approaches dis- 
cussed in previous sections, This section shows results obtained from a simulation 
model. The goals of the presentation is to illustrate the theoretical results. Finally, 
a flight teat data set is studied in the next section to qualitatively evaluate the need 
for indicial terms. 


5.i Simul atio n _D&La 

The simulated model 13 of the form 


q ~ M„a + M.,q -I- Mt,,b c -I- j Mq{r)q[t - r)cir (32) 

0 

A dimensional model has been selected to eliminate various scaling parameters. The 
following values are used in simulations, throughout 

Zo = -2 , M„ = -5 , A/, = -i , M t , * —5 (33) 

Note that we have selected an indicial term in q , a linear combination of a and a. 
Of course, this term can also be written in terms of a, using integration by parts. 
The characteristic equation of this mode' ‘s j 2 -Ms + 9 . So the natural frequency is 
3 rad sec (about | Hz), 

In our simulations, wc chose two values of T, .5 and .0“ , uiT then is 1.5 and .15. 


ujT :.s 

The modeled form of M^(t) is shown in Figure 2, The simulation model was 
developed rapidly using SystemJ3uild , a feature of MATRIX x (10). The block 
diagram is shown in Figure 3. An input time history shown in Figure 4 is applied 
to give an a , q , q response of Figure b. 

At first least squares estimation approach is used and < j is derived by diffeienti- 
at'.ng q (a simple two point differentiating is used throughout). The . esulting model 
without noise is shown in Figure 6. Almost all of the error is due to uncertainties in 
obtaining <? from q . The estimated parameter values are exact if simulation values 
of (j is used 1 :ect!y in estimation. 

Figure 7 shows the effect of adding noise to q on model estimation accuracy. 
Table 1 snows theoretical values of estimation error for random noise in q . These 
results indicate that the least squares approach does not give an acceptable estimate 
for indicial models unless a more complex input is used. The least squares method 
nevertheless is good for selecting significant terms. 

The approximation to this indicial model would be as follows. 

7 - c a- 2q~ 56' - 2.1, (34) 
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When 'an attempt is made to identify tliis ; model,. 'tlie estimated parameters are 
shown in- Table 2 and the fit error to ? is shown 1 id Figure 8. ? 

Maximum likelihood technique is used to estimate the indicial model from noisy 
measurements of angle-of-attack, a, and pitch rate q. The estimated indicial function 
is shown in Figure 9, and the fit to a and q time history is shown in Figure ,10. 

wr = 0.15 ■ " 

The indicial function is now spread out over the first .05 sec, but is scaled up by 
a factor of 10 from Figure 2, such that the approximation is the same as in Equation 
(34). The input is the same as before. The responses are shown in Figure 11. The 
following steps are repeated. 

(1) q is obtained by differentiating q and the indicial model is identified. The 
results are shown in Figure 12 and are much poorer than before. 

(2) The estimates are so poor without noise that very little noise causes the 
estimates to lose all accuracy. 

(3) Table 3 shows that the simplified model may be identified better than for 
w T = 1.5 Figure 13 shows fit errors in q . 

(4) Application of maximum likelihood approach does not aid the estimation of 
the indicial function. 

5.2 Flight Data 

This section shows data from a spin research vehicle (SRV) and demonstrates 
the need to use indicial function representations in the high angle-of-attack region. 
The 3/8-scaled unpowered model of a high performance fighter aircraft is dropped 
from a B-52 and is controlled remotely. The vehicle provides an effective means to 
conduct aerodynamic tests in post-stall and spin regions. Though the flight test 
lasts several minutes, a 40-second long segment shown in Figure 14 is studied. A 
standard set of on-board instrumentation is available including (3-axis rate and 
attitude gyros and 3-axis accelerometer, dynamic pressure, angle-of-attack, sideslip 
angle, static pressure, and control position). 

We will look at data qualitatively in this report to show the need for integro- 
differential function models or a function with memory. Figure 15 shows the normal 
force coefficient time history as a function of angle-of- attack. Based on known air- 
craft behavior, the lift coeffi ’ent depends primarily on angle-of-attack and to a 
lesser extent on elevator deflection and other aerodynamic variables. The depen- 
dence of Cw on variables other than the angle-of-attack cannot account for the large 
loops. Note that for fixed controls between 28s and 40s, the aircraft exhibits os- 
cillatory behavior with increasing angle-of- attack. This is an interesting region to 
analyze because forces and moments are affected only by aerodynamic variables. 

Figure 15 indicates that the phenomenon i3 not aerodynamic hysterisis, because 
in hysterisis the function follows one path for increasing and another one for de- 
creasing values of the independent variable. The behavior of the aircraft for angle- 
of-attack between 70° and 80° is qualitatively very interesting. It is known from 
linearized analysis [8] that over 40 6 angle-of-attack, the magnitude of C Ntt decreases 
with angle-of-attack. In the time history shown in Figure 16 the angle-of-attack 
suddenly jumps from about 40° to 80" and the lift stays high (corresponding more 
to the 40° lift coefficient than the 80° lift coefficient). Then about 0.5 sec later the 
lift coefficient drops suddenly. PhyscialJy, this may occur because a suddon increase 
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in anglc-of-attack causes the flow to remain attached for a short period of time and 
the lift increases dramatically. The flow must events. -lly separate causing the lift to 
drop. The phenomenon may be related to what has been referred to as “dynamic 
lift". 


The plots also show that the natural frequency is about 0.5 Hz. Thus an indicia! 
function with 0.5 sec. delay should be clearly identifiable. The identified model had 
been shown previously as Figure 17 in Reference [ 5] . 

Note that for responses below 40° angle-of-attack, the indicia! model is not 
identifiable and the indicial term may be approximated by a lumped coefficient. 
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VI. SUMMARY AND CONCLUSIONS 

The report presented approaches to identify intergro-different.al equation mod- 
els which arise in aeroelastic and unsteady aerodynamic.:. The following conclusions 
may be drawn. 

1 . When the product of frequency of motion and maximum time delay is much 
smaller than one, the integral term can be approximated by a constant. 
When this product is of the order of, or larger than one, the integral term 
cannot be approximated, 

2. lntegro-differential models are in general non-identifiabie. Approximations 
are needed to bring about iclentifi ability. 

3. Least-squares method may be used foi model determination but the maxi- 
mum likelihood technique is needed to accurately estimate parameters. 

4. High angle-of-attack and post-st.all/spin region appears to have character- 
istics, which can be satisfied by indicial models. 

More work though is needed to advance better understanding of unsteady- aerody- 
namics and aeroelastic phenomena from measured data. 
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Appendix A 


APPROXIMATION OF J ff mA (r)a(t - r)dr 


o 


We will attempt to approximate this integral by C nla , For this approxiamtion 
to hold in steady state 


<?,„„(() = // r „ 4 (oo) 

The error is approximated as follows: 


(A.l) 


c ~ J H m A [r)a(t - r)dr - C mo a(t) 

o 

CO 

- J {H m/ . (’■) - ffma (oo)}a(t - r)dr (4.2) 

o 


It is reasonable to assume that //„ lA (r) reaches its steady-state value // mA ( oo) for 
r T , Then 


|//m A (r) - U mA (oo)| < (T ~ r)\H mA \ mtlx r < T 
« 0 r > T 


(A3) 


Then (A. 2) becomes 


e<c|//„J max &tnax T 2 (A.4) 

where c is a constant. Since H m<i is the same as F ma (r) and the integral of F rrl „(r) 
from 0 to T is Equation (A. 4) becomes 


£ ^ c ot max T 


(- 4 . 5 ) 


Equations can also be developed to approximate the integral term by a combi- 
nation of a and a coefficients. 
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TABLE L: COMPUTED ESTIMATION ERROR FOR LEAST SQUARES 
(500 Points, RMS Noise in q = .001) 


Parameter 

Value 

Estimation Error 

M„ 

-5.0 

.0048 

M,, 

-2.0 

.0167 

M„. 

-5.0 

.0045 

M„ 

-0.0 

.00004 

M,(. 1) 

-0.3 

.0253 

AM- 2 ) 

-0.6 

.0213 

M„(- 3) 

-0.7 

.0183 

AM-4) 

-0.4 

.0128 

AM-5) 

-0.2 

.0045 


TABLE 2: ESTIMATED PARAMETERS IN REDUCED MODEL 
(wr = l.g) 


Parameters 

Values 

Estimated Value 
(No noise) 


-5.0 

-9.75 

M,,[equi.) 

-4.1 

-3.79 


-5.0 

-6.28 


0.0 

fO.O 
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TABLE 3: ESTIMATED PARAMETERS IN REDUCED MODEL 

(wT .15) 


Parameters 

Values 

Estimated Value 
(No noise) 

Ma 

-5.0 

-5.3 

M,,(equx.) 

-4.1 

-4 A 

Mt. 

-5.0 

-5.3 

M 0 

0.0 

- 1 - 0.0 
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Figure 1 
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: Schematic Procedure for the Estimation of Integro- Differential Models. 



Figure 2 : Indicial Models Used in Simulation 
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Figure 7 : Effect of Noise in q on Least Squares Estimation Accuracy 
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Figure 8 : Fit-Error in g when an Indicial Model is not Used in Estimation (uT - 1.5) 



























Figure 15 


Figure 16 : 
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: Normal Force Coefficient as a Function of Angle-of-Attack 
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Figure 17 : Estimated Indicial Function 
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